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We present a measurement of the time-dependent CP asymmetry for the neutral B-meson decay 
B° cj>K°. We use a sample of approximately 114 million B-meson pairs taken at the T(4S) 
resonance with the BABAR detector at the PEP-II B-meson Factory at SLAC. We reconstruct 
the CP eigenstates (f>Ks and 4>Kl where (f) K~ , Ks tv'^tt', and K1 is observed via its 

hadronic interactions. The other B meson in the event is tagged as either a B" or B^ from its decay 
products. The values of the CP-violation parameters are S^K = 0.47±0.34(stat)+oog(syst) and 
C^K = 0.01 ± 0.33(stat) ± O.lO(syst). 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



Decays of B mesons into charmless hadronic final 
states with a (/> meson are dominated by 6 — *■ sss gluonic 
penguin amplitudes, possibly with smaller contributions 
from electroweak penguins, while other Standard Model 
(SM) amplitudes are strongly suppressed 0. In the 
SM, CP violation arises from a single complex phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing 
matrix 0. Neglecting CKM-suppressed contributions, 
the time-dependent CP-violating asymmetries in the de- 
cays — > (j>K° and P" J/ipK° are proportional 
to the same parameter sin 2/3 where the latter de- 
cay is dominated by tree diagrams. Since many scenar- 
ios of physics beyond the SM introduce additional dia- 
grams with heavy particles in the penguin loops and new 
CP- violating phases, comparison of CP- violating observ- 
ables with SM expectations is a sensitive probe for new 
physics. Measurements of sin2/3 in B decays to charmo- 
nium such as P" — > JhbK^ have been reported by the 
BABAR ^ and Belle collaborations, and the world 
average for sin2/3 is 0.731 ± 0.056 The Belle collabo- 
ration measures sin2/3 = —0.96 ± 0.501;q'5i in the decay 
B"^(j)K°^. 

In this letter we report a measurement of the time- 
dependent CP asymmetry in the final state ° based 
on an integrated luminosity of approximately 108 fb"^ 
collected at the T(4S') resonance with the BABAR detec- 
tor at the PEP-II asymmetric e+e" collider located 
at the Stanford Linear Accelerator Center. 

From a B^B^ meson pair we fully reconstruct one me- 
son, Bcp, in the final state <j}K'^ , and partially recon- 
struct the recoil B meson, Ptag- We examine Ptag for 
evidence that it decayed either as P° or B^ (fiavor tag). 
The asymmetric beam configuration in the laboratory 
frame provides a boost of (3^ = 0.56 to the T(4S'), which 
allows the determination of the proper decay time differ- 
ence At = tcp — ^tag from the vertex separation of the 
two neutral B mesons along the beam [z) axis. The de- 
cay rate f+(f-) when the tagging meson is a B^{B'^) is 
given by 

f±(Ai) = — [l± S^Ksin{AmdM) 

T C^K COS {Anid At)], (1) 

where Tgo is the neutral P meson mean lifetime, and 
Arrid is the B'^-B^ oscillation frequency. The time- 
dependent CP- violating asymmetry is defined as Acp = 



(f+ — f_)/(f+ + f-). In the SM, decays that proceed 
purely via the h —f sss penguin transitions have CP 
parameters S'^if = —r// sin 2/3 and C^k = 0, where 
/3 = engi-V^^V^JV^^V*]. Here Vik is the CKM matrix 
element for quarks i and fc, and the CP eigenvalue is 
rjf = -1 (+1) for </.i^o i^Kl). 

The Be p candidate is reconstructed in the decay mode 
(/iX° with (t) the K° is either a K° or a 

if 3 — > 7r+7r~. We combine pairs of oppositely charged 
tracks extrapolated to a common vertex to form (/) and 
K'^ candidates. For the charged tracks from the decay 
we require at least 12 measured drift-chamber (DCH) 
coordinates and a minimal transverse momentum px of 
0.1 GeV/c. The tracks must also originate within 1.5 cm 
in xy and ±10 cm along the z-axis of the nominal beam 
spot. Tracks used to reconstruct the </> mesons are distin- 
guished from pion and proton tracks via a requirement 
on a likelihood ratio that combines dE/dx information 
from the silicon vertex tracker (SVT) and the DCH for 
tracks with momentum p < 0.7GeV/c. For tracks with 
higher p, dE/dx in the DCH and the Cherenkov angle and 
the number of photons as measured by the internally re- 
fiecting ring-imaging Cherenkov detector are used in the 
likelihood. The two-kaon invariant mass must be within 
16 MeV/c^ of the nominal cj) mass 6]. 

For tracks corresponding to if^ and Ptag daughters our 
requirements are less restrictive. A 7r~''7r~ candi- 

date is accepted if its two-pion invariant mass is within 
15 MeV/c^ of the nominal mass 0, its reconstructed 
decay vertex is separated from the collision point by at 
least 3 standard deviations, and the angle between the 
line connecting the (j) and decay vertices and the 
momentum direction is less than 45 mrad. 

We identify a candidate as in our P° J/ipK^ 
analysis either as a cluster of energy deposits in the 
electromagnetic calorimeter (EMC) or as a cluster of hits 
in two or more layers of the instrumented fiux return 
(IFR) that cannot be associated with any charged track 
in the event. The energy is not well measured. There- 
fore, we determine the laboratory momentum from 
its fiight direction as measured from the EMC or IFR 
cluster and the constraint that the invariant (pK^ mass 
agrees with the known P° mass. In those cases where 
the K° is detected in both the IFR and EMC we use the 
angular information from the EMC, as it has a higher 
precision. In order to reduce background from tt*^ de- 
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cays, we reject an EMC candidate cluster if it forms 
an invariant mass between 100 and 150 MeV/c^ with any 
other cluster in the event under the 77 hypothesis, or if it 
has energy greater than 1 GeV and contains two shower 
maxima consistent with two photons from a 7r° decay. 
The remaining background of photons and overlapping 
showers is further reduced with the use of a neural net- 
work constructed from cluster shape variables, trained 
on Monte Carlo (MC) simulated B° 4)Kl and mea- 
sured radiative Bhabha events, and tested on measured 
e+e" ^ <l){-^ K°K°)j and B° J/^K^ events. The fi- 
nal (t)Kl sample consists of approximately equal numbers 
of IFR and EMC candidates. 

The results are extracted from an extended unbinned 
maximum likelihood fit for which we parameterize the 
distributions of several kinematic and topological vari- 
ables for signal and background events in terms of prob- 
ability density functions (PDFs) The background 
arises primarily from random combinations of tracks pro- 
duced in events of the type e+e~ qq, where q — 
u,d,s,c (continuum). Background from other B decay 
final states with and without charm is estimated with 
MC simulations. Opposite CP contributions from the 
K~^K~K'^ final state (e.g. K^K'^ S-wave) are estimated 
with an angular analysis on data to be less than 6.6% and 
treated as systematic error. The shapes of event variable 
distributions are obtained from signal and background 
MC samples and high statistics data control samples. 

Each Bcp candidate is characterized by the en- 
ergy difference A£' = Eg — and, except for 
B° cl)Kl, the beam -energy-substituted mass mEs = 

\J {\3 +pa ■ pbY /El — p\ "s]. The subscripts and B 

refer to the initial T(4S') and the Bcp candidate, respec- 
tively, and the asterisk denotes the T(4S') rest frame. For 
signal events, AE is expected to peak at zero and tties 
at the nominal B mass. We require /S.E < 0.08 GeV for 
B° (l)K° and |A£;| < 0.2 GeV and uies > 5.2GeV/c^ 
for (pKs- the fit we also use the helicity angle 

9h, which is defined as the angle between the directions 
of the and the parent Bcp in the K^K~ rest frame. 
The cos Oh distribution for pseudoscalar- vector B decay 
modes is cos^ Oh, and for the combinatorial background 
it is nearly uniform. 

In continuum events, particles appear bundled into 
jets. This topology can be characterized with several 
variables computed in the CM frame. One such quantity 
is the angle 9t between the thrust axis of the Bcp candi- 
date and the thrust axis formed from the other charged 
and neutral particles in the event. We also use the an- 
gle 9b between the Bcp momentum and the beam axis, 
and the sum of the momenta pi of the other charged and 
neutral particles in the event weighted by the Legendre 
polynomials Ln{9i),n = 0,2, where 9i is the angle be- 
tween the momentum of particle i and the thrust axis 
of the Bcp candidate. For i?'^ candidates, we 



combine these variables into a Fisher discriminant [l^l 
after requiring |cos0t| < 0.9. In this mode background 
from other B decays is negligible, as demonstrated in MC 
simulation studies. 

More stringent criteria must be applied to supress 
backgrounds in the case of B'^ 4>Kl candidates, and 
we require |cos0t| < 0.8 and |cos0b| < 0.85. We de- 
fine the missing momentum Pmiss , calculated in the lab- 
oratory frame from the beam momentum and all tracks 
and EMC clusters, excluding the candidate. We re- 
quire the polar angle 9miss of the missing momentum 
with respect to the beam direction to be greater than 
0.3 rad. The cosine of the angle between Pmiss and the 
direction, 9k, must satisfy cos9k > 0.6. In the plane 
transverse to the beam direction, the difference between 
the missing momentum projected along the direc- 
tion and the calculated momentum must be greater 
than —0.75 GeV/c. In the Fisher discriminant we replace 
I cos 05 1 by the cosine of the angle between the missing 
momentum and the from the (p decay. In the (j>Kl 
sample about 1.4% of the events originate from charm 
B decays, 0.7% originate from charmless B decays, and 
about 0.2% potentially have a CP asymmetry. The dom- 
inant contamination is the mode B — s- (pK* , where the 
K* decays to K'^tt. In the likelihood fit we explicitly pa- 
rameterize backgrounds from both charm and charmless 
B decays as derived from MC simulations. 

All the other tracks and clusters in the event are used 
to form the -Btag, and its flavor is determined with a mul- 
tivariate tagging algorithm The tagging efficiency e 
and mistag probability w in four hierarchical and mu- 
tually exclusive categories is measured from fully recon- 
structed B° decays into the D^*'>-X+ {X+ = TT+,p+,af) 
and J/%l}K*^{K*° K+tt-) flavor eigenstates (Bflav 
sample). The analyzing power e{l — 2w)^ is (28.7±0.7)%. 

A detailed description of the At reconstruction algo- 
rithm is given in Ref. The Bcp vertex resolution 
is dominated by the (j> vertex. The average Az resolu- 
tion is 190 fim. and is dominated by the tagging vertex 
in the event. Thus, we can characterize the resolution 
with the much larger Bflav sample, which we fit simul- 
taneously with the CP samples. The amplitudes for the 
Bcp asymmetries and for the Bflav flavor oscillations are 
reduced by the same factor due to wrong tags. Both 
distributions are convoluted with a common At resolu- 
tion function, and the backgrounds are accounted for by 
adding terms to the likelihood, incorporated with differ- 
ent assuinptions about their At evolution and resolution 
function [lOj. 

Since we measure the correlations among the observ- 
ables to be small in the data samples entering the fit, 
we take the probability density function for each 
event j to be a product of the PDFs for the sepa- 
rate observables. For each event hypothesis i (signal, 
background) and tagging category c, we define Vf^ = 

(mEs) -V^ (AE) -r, {T) ■V^ (cos 9H)■V^ {At; a At , c) , where 
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for the (j)Kl mode Vi{mEs) = 1 and for the flavor sample 
'Pi{T) ■ Vi{cos9H) = 1- The a At is the error on At for a 
given event. The likelihood function for each decay chain 
is then 



E 



(2) 



where Ni^c is the yield of events of hypothesis i found by 
the fitter in category c, and Nc is the number of category 
c events in the sample. The total sample consists of 86200 
Bflav, 2138 (t>K° and 4730 (l)Kl candidates. We find 70±9 
(j)Kg and 52 ± 16 (j)Kl signal events. The signal yields in 
both the (pK'^ channels agree well with our determination 
of the branching fraction for B° (j}K° [l^. Fig. [2 
shows the tties {AE) distribution of (j)K'^ {(f'K^) events 
together with the result from the fit after a requirement 
on the likelihood (computed without the variable plotted) 
to enhance the sensitivity. 




5.22 5.24 5.26 5.28 5.3 
mEs[GeV/c^] 



FIG. 1: Distribution of the event variable (a) niES for the 
4>Kg final state and (b) AE for the (j)Kl final state after 
reconstruction and a requirement for the likelihood with total 
signal efficiency of 32% and 5%, respectively. The solid line 
represents the fit result for the total event yield and the dotted 
line for the background. 



We determine the CP parameters S^k and C^k along 
with an additional 38 free parameters: the efficiency per 
tagging category (4 parameters), the average mistag frac- 
tion and the difference between i?*^ and mistags for 
each tagging category (8 parameters), the signal At res- 
olution (9), and time dependence (6), At resolution (3) 
and mistag fractions (8) for the background. We fix 
Tbo and Am^ to the world averages The determi- 
nation of the mistag fractions and Ai-resolution param- 
eters is dominated by the high-statistics iSflav sample. 
The fit was tested with a parameterized simulation of 
a large number of data-sized experiments and full de- 
tector simulated events for the different signal and back- 
ground samples. The fit was also verified with our J /tpKl 
data sample and a control sample of 232 (f)K'^ candi- 
dates where one expects 5^^+ = C^k+ = 0. We mea- 
sure S^K+ = 0.23 ± 0.24 and C^k+ = -0.14 ± 0.18 with 
statistical errors only. The simultaneous fit to the 4)K^ 




FIG.^: Plots a) and b) show the At distributions of _B°- 
and _B"-tagged events. The solid lines refer to the fit 

for all events; the dashed lines correspond to the background. 
Plot c) shows the asymmetry. A requirement for the event 
likelihood is applied. Plots d), e), and f) are the corresponding 
plots for 4>K1 events. 



and flavor decay modes yields: 

S^K = 0.47±0.34(stat)t°:°^(syst), 
C^K = 0.01±0.33(stat)±0.10(syst). 

The result in the dominant channel _B" — > 4>K^ is S^k = 
0.45±0.43 and C^if = — 0.38±0.37 with statistical errors 
only. Fig.|21shows the At distributions of the B^- and the 
-B^-tagged subsets together with the raw asymmetry for 
the (i>K^ and (t>K^ events with the result of the combined 
time-dependent CP-asymmetry fit superimposed. 

We consider systematic uncertainties in the CP coeffi- 
cients S(j,K and Ccj,K due to the event-yield determination 
in the two channels (±0.01 for S^k, ±0.05 for C^k), con- 
tributions from B^ final states with opposite CP (+0.06, 
±0.02), the parameterization of PDFs for the event yield 
in signal and background (±0.02, ±0.05), composition 
and CP asymmetry of the background in the CP events 
(±0.03, ±0.03), the assumed parameterization of the At 
resolution function (±0.02, ±0.01), the mss background 
parameterization (±0.02, ±0.05), a possible difference in 
the efficiency for B° and B° (±0.01, ±0.02), the fixed 
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values for Am^ and tb (±0.00, ±0.01), the beam-spot 
position (±0.01, ±0.01), and uncertainties in the SVT 
alignment (±0.01, ±0.01). The bias in the coefficients 
due to the fit procedure (±0.03, ±0.01) is included in 
the uncertainty without making corrections to the final 
results. We estimate errors due to the effect of doubly 
CKM-suppressed decays Q to be (±0.01, ±0.03). We 
add these contributions in quadrature to obtain the total 
systematic uncertainty. 

In summary, we have measured the time-dependent CP 
asymmetries in the combined i?-meson final states 4*^3 
and (pK^. We obtain values for the CP- violation param- 
eters S^K and C^if that agree within one standard devi- 
ation with the ones measured in the charmonium chan- 
nels 0,01; the central value of S^k is also consistent with 
no CP asymmetry at the 1.3 cr level. 
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